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Nomenclature 


English letters 

a = indicates air 
c = chord length 
D = drop 

f = subscript, denotes fog or vapor state 
i = subscript, denotes incident 
L = characteristic length 
t = time 

U, = free stream velocity 

u,v = velocity components along x-y-directions respectively 
w = subscript, denotes water 
x,y = physical coordinates 


Greek letters 

a = angle of attack, or non-dimensional ration = W ,/p 
P = droplet incident angle ^ 

p = density of the fluid 
Pv = density of the vapor 
pi = density of the liquid 
p = coefficient of viscosity 
5 = thickness of the boundary layers 
v = kinematic coefficient of viscosity = p/p 
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SUMMARY 


A mathematical model for a two-phase flow laminar airfoil in simulated heavy 
rain has been established. 

The set of rion-linear partial differential equations has been converted into a set 
of finite difference equations; appropriate initial and boundary conditions are 
provided. The numerical results are compared with the experimental 
measurements. They show good agreement in quality. 


I. INTRODUCTION 

This paper investigates the rain effects to the airfoil during landing and take-off. 
When the heavy rain hits the Airfoil, the fine water drops near to the wing surface 
form a liquid film while the water vapor above the wing surface establishes a 
gaseous fog. The former is in liquid state, the latter is in gaseous state. This is the so- 
called two-phase phenomenon. The thin liquid film close to the wing surface 
usually forms a laminar boundary layer. However, the heavy rain fall may cause a 
premature turbulent boundary layer. That is, the heavy rain effects may cause the 
laminar boundary layer transit into turbulent boundary layer. As a result, the lift is 
appreciably decreased, while the drag is considerably increased. In a a more serious 
situation, the decrease in lift and increase in drag may develop into wind shear 
which causes an airplane catastrophe involving loss of human life and property 
damage. This important knowledge should keep the pilot well-informed as a 
precaution of flight safety. A variety of experimental and analytic methods used to 
investigate the heavy rain effects on airfoil may be referred to in previous works, 
notably references # 1 through # 13. This paper has established a mathematical 
model for the airfoil under the simulated heavy rain. The said mathematical model 
consists of a set of non-linear partial differential equations for which a numerical 
solution has been obtained. A comparison between theory and experiment has been 
made and will shed some light on flight safety. 

n. FORMULATION 


As in the following figure, the two-phase-flow around an airfoil is clearly 

RAINDROPS wteractwq with an airfoil 


FIGURE (1) 
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described. To simplify the problem the following assumptions are made: 

1. The fluid flow is non-steady, viscous, and incompressible. 

2. There is laminar boundary in the flow region. 

3. The airfoil is represented by a flat plate. The physical coordinates are shown 
in Figure (2). 

For a two-dimensional, nonsteady two-phase flow, we obtain the following 
sets of fundamental equations under boundary layer approximation: For the liquid 
phase 


Equation of Continuity: 


Equation of 2 u + 

Momentum : ^ C 
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(2) 


Film thickness: 5 = 5-^ux /U 2 


For the fog phase Equation of Continuity: 


Equation of Momentum: 


du 

dx 






(3) 

9u + dv 
dx dy 

= 0 

... ... 

(4) 

du 9u 

J. XT’ _ 

= M_ 

a 2 u 
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dx dy 

Pr 
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The initial conditions are: 


at t = 0, u = U (x,0), V = 0, 5 = 0 
2 

The boundary conditions are: 

y = 0, for the liquid phase 
u = 0, v = 0 (no slip condition) 

at y = 5 (liquid-fog interface) 

(u)l = (u)f (v)l = (v)f 


( 3u^ i 

3u'| 

II 
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In other words, U, V, and 3u/3y for both liquid and fog must be compatible at the 
interface, as y -»«> (for the fog phase ) t u - U(x,t) (9) 


m. NON-DIMENSIONALIZATION 


In order to non-dimensionalize the set of fundamental equations, we 
introduce the following: 


x 

5 


x 


L 

_5 

L 



u 






( 10 ) 


where Uc is the reference velocity, and L is the chord length of the airfoil. 
Substituting the above relations into Equations (1), (2), (3), (4) and (5), and omitting 
the bar notations, we obtain: 
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( 11 ) 

( 12 ) 

(13) 


In Equ. (12), Wi has the dimension of M/L^, so does the density p. Wi/pw is a 

non-dimensional quantity, where (RN)i = p w U2x/ p. Reynolds number of the liquid 

water. We let a = Wj/ p w = non-dimensional ratio of the liquid water content to the 
water density. 

For the fog phase 


Equation of Continuity: 
3u 3v 

37 + 37 = 0 - ■ 

Equation of momentum 


( 14 ) 
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(15) 


9u du 5u 1 9 2 u 

— + u h v- — = — T 

dx dx dy (RN) f 3y 2 


where (RN) ( = = Reynolds number of the fog. 

The corresponding initial and boundary conditions are: 
atx = 0, u=U(x,o), V = 0 

5 = 0 



(6-a) 


for tiie liquid phase, at y = 0, u = v=o (No slip condition) ... .... (7.a) 

at y = 5/L (Liquid-vapor interface) 


H = 

(u), (V), = (V), 


f 3u l 

f du) 

1 (8. a) 


= HT- 

l d y J ( 


As y — » « 

u = U(x,t) ... ... 

(9 a) 


IV. FINITE DIFFERENCE EQUATIONS 

The previously derived set of non-dimensional partial differential equations can 
be transformed into the finite difference equations in the following manner: An 

explicit method is used. Let U’, V' and 6’ denote the values of U, V and $ at the 
end of a time-step. Then the appropriate sets of finite difference equations are: 

For the liquid 

Equation of Continuity: 
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(ll.a) 


<1 - <u , V» ~ V?-m 

Ax Ay 

Equation of Momentum: 

U M - P. |TJ ( ^ ) 

Ax “ Ax 
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1 
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Film thickness: 


5 U - 5 j 


vIAx 

U 2 
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For the fog 

Equation of Continuity: 


K - Hk> 

Ax 


v m - ZLli 

Ay 


= 0 


Equation of Momentum: 

tU ( u u - u,., f ) 

At Ui '* Ax 


+ V, 




Ay 


_ 1 (U w - zu u + U,. M ) 

(RN), 

The boundary conditions are 

at x = 0, u = U 2 , v = 0 I 
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(12. a) 


(13.a) 


(14.a) 


(15. a) 


(16) 


for the liquid phase, y = 0, U=V=0 .. (17) .At the interface 0<x<L,O< y 


5max 
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(16.a) 


(16-b) 


All the points denoted by Equations (16.a) and (16.b) are boundary points at which 
the values of U and V are already known. Furthermore the shearing force of the 
interface is denoted by the following relation 


(U. . , u ) 

V1,J+1 - i.j 


\(L * V/ 

A y 


a y 


(17.a) 


V. THE STABILITY OF THE FINITE -DIFFERENCE EQUATIONS 


Since an explicit procedure is used, we wish to know the largest time-step 

consistent with stability. Equation of continuity is ignored since Ax does not appear 
in it. The general terms of the Fourier expansion for U at a time arbitrarily called 

x=0 are both e^ ax e^Py, apart from a constant (Here i = V-l). At a time x later, these 
terms will become 


U: \|f(x) e iax *e i(Jy 

Substituting the above into Equation (12.a), regarding the coefficients U and V as 
constants over any one step, and denoting the values after time-step by y' gives 

2j2t ♦ + y jjgm 

Ax Ax Ay 


l 2vji(x)[cos((3Ay)-l] a /2 

RN ( A y ) 2 


+ -T- V, 2 sinpcosPe- i(ax+(Jy) 
5, > 




(18) 
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Through a very tedious algebraic manipulation, there obtains the criterion of 
stability: 


U y +_2Ax_ 

Ax Ay RN(Ay) 2 


£ 1 


(19) 


In the present research, Rn ,the Reynolds number, is in the order of ICR 
Equation (19) follows automatically. However, the coefficients U and V, treated as 
constants over any one time-step, will vary from one time-step to the next in a 
manner which cannot be predicted as a priori. That is the maximum permissible 
time-step consistent with stability and is itself variable, but its value can always be 
checked during computation if necessary. 


VI. CALCULATION OF LIFT COEFFICIENT 
According to Glauert, the lift coefficient is given by 


CL = rc(a+£ 0 ) (20) 


where d is the angle of attack. 



( 21 ) 


and f (x/c) = 
1 



The relationship between (x/c) and f 


(x/c) is : 


( 22 ) 


(x/c) 

0.025 

0.05 

0.10 

0.20 

0.30 

0.40 

f 1 (x/c) 

2.090 

1.54 

1.18 

1.00 

0.99 

1.08 

(x/c) 

0.50 

0.60 

0.70 

0.80 

0.90 

0.95 

fi (x/c) 

1.27 

1.62 

2.31 

3.98 

10.60 

29.20 


For Wortmann airfoil, the value of (y/c) for each (x/c) can be found. The 
integration can be performed easily. However, for flat plate airfoil, £o = O, the lift 
coefficent reduces to Cl = xa. 
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yn. CALCULATION OF DRAG COEFFICIENT 

A body of buff form, particularly if it has sharp edges like a flat plate (which is 
the case here) will shed strong vortices and will have a large form drag. As a first 
approximation, the form drag coefficient may be written 


C D 



•— = 0.281 t~r 

V bV 


(23) 


where K = I'll a U/ is the strength of each point vertex, the so- called Karman vortex 
street, a is the distance separating the successive vortices of each row, U is the 
induced velocity, h is the distance between the two rows, i.e. the breadth of the 
street. A sketch for von Karman vortex street is given as follows : 
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Figure (3). 
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The skin frictional drag coefficient of a flat plate is given by a 


classical formula: 


Cd- 1 


.zizjfW 


The historical results and the current computation are tabulated as follows : 



As shown in the above table, the drag coefficient of the airfoil with 
is higher than that without rain. The computed results are much 
lower than those measured by Hansman(in references 5, and 6). The reason for 
this is that the over-simplified mathematical model of an airfoil by a flat 
plate does not cover the physical reality of a true airfoil which has proper 
thickness and suitable chord length. 

Also the premature transition of laminar flow into turbulent flow may 
cause higher drag coefficient . This possible turbulent phenomenon is not 
covered due to our previous assumptions. 

The superposition of form drag to the skin frictional drag should improve 
the agreement in quantity with the experimental measurement. However, the 
analytic determination of form drag in Eq. (23) remains to be a difficult 
task for further investigation. 
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VIII Condusion and Discussion 


The set of finite difference equations has been converted into Fortran language, 
and a numerical solution is obtained. The complexity and computation time are far 
beyond the investigator's antidpation. For the numerical results it was found that 
there is decrease in lift and increase in drag due to the heavy rain effects. The reason 
for this is that the heavy rain causes roughness on the wet surface of the airfoil. A 
comparison of the present numerical calculation with Hansman's (6) experimental 
measurement is shown in Figure (5). According to Hansman, the Wortmann 
section had the greatest lift degradation: nearly 25%. The computed decrease in lift 
and increase in drag are both lower than that of the experimental measurements. 
The reason may be that the over-simplified mathematical model could not cover 
the physical reality, such as the premature boundary layer transition near the 
leading edge of the airfoil and the three-dimensional effects. 

The laminar layer thickness of both the liquid film and fog is shown in Figure 
(6). The velodty profile near the surface of the airfoil is shown in Figure (4). It is 
evident that the U-component of the velocity near the surface of the airfoil is 
decreased due to the rain-effects. 

It requires further investigation for a turbulence model and three-dimensional 
wing in order to accomplish perfect agreement between analytical investigation and 
experimental measurements. 
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- velocity profile near the surface of an airfoil interacting with rain 
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Angle of Attack 


Figure (5). Lift and drag coefficients vs angle 

of attack for the Wortmann FX67-K170 
airfoil in dry and wet conditions. 
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